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Abstract. This article studies multiple scattering of matter waves by a disordered optical 
potential in two and in three dimensions. We calculate fundamental transport quantities such 
as the scattering mean free path the Boltzmann transport mean free path and the 
Boltzmann diffusion constant Db, using a diagrammatic Green functions approach in the 
weak scattering regime. Coherent multiple scattering induces interference corrections known 
as weak localization which entail a reduced diffusion constant. We derive the con'esponding 
expressions for matter wave transport in an correlated speckle potential and provide the 
relevant parameter values for a possible experimental study of this coherent transport regime, 
rncludrng the critical crossover to the regime of strong or Anderson locaUzation. 



1. Introduction 

The observation of the first gaseous Bose-Einstein condensates in 1995 and, a few years later, 
of the first gaseous ultra-cold fermion gases constitutes a major achievement in the field of 
atomic physics. Loading ultra-cold degenerate gases, be it fermions, bosons, or fermion- 
boson mixtures into optical lattices has opened fascinating new perspectives for the study 
of condensed-matter quantum physics [1^]. A particularly promising line of research is 
the experimental and theoretical study of quantum phase transitions [5]. Indeed, the high 
degree of control and precision achieved in these experiments has allowed systematic studies 
of physical phenomena that are observed until now mostly in solid-state systems. 

A natural evolution in this domain is the investigation of the influence of disorder [6] 
which can induce a Bose-glass phase for strongly interacting bosons [7], and a Lifshitz 
glass [8] or an Anderson glass [9] for weaker interaction. Experimentally, a major milestone 
has been reached with the realization of quasi- ID condensates evolving in a speckle light 
field [10-14]. Such disordered optical potentials can be easily generated and their statistical 
properties are well known [15]. In two dimensions, atomic diffusion has been studied in 
optical quasi-crystals with five-fold symmetry for atomic clouds in the dissipative regime [17] 
and Bose-Einstein condensates [18] covering the intermediate regime between ordered and 
completely disordered systems. Most recently the experimental observation of the onset of the 
Bose glass phase for ultra-cold atoms in a bichromatic optical lattice has been reported [16]. 
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Disorder has long proven to be a crucial ingredient to understand coherent transport 
properties. The most prominent example is the weak localization phenomenon in mesoscopic 
physics [19-21], which has been studied extensively for electrons [22,23] and for classical 
waves [24, 25]. Weak localization arises from interference between multiply scattered waves 
in a random medium. This interference survives the configuration average over many 
realizations of disorder and reduces the conductivity that enters the Drude model for electron 
transport [26] and the diffusion constant of classical radiative transport theory [27, 28]. 
Disorder can even induce a metal-insulator transition, known as the strong (or Anderson) 
localization transition [29-33], that has been studied extensively in the framework of the tight 
binding model [34, 35]. In the localized regime the quantum states cease to extend over the 
whole system and become spatially localized with the consequence that quantum transport 
through the system is exponentially suppressed. 

From a theoretical point of view, atom transport shares many similarities with the 
radiative transfer theory [36] and even more with electron transport theory [20, 21]. 
Unfortunately, the unambiguous observation of wave localization phenomena for matter 
waves is difficult if the quantum evolution is interaction-dominated. But there are several 
ways to achieve the interaction-free regime: for example, one can let the ultra-cold gas 
first expand, so as to decrease its spatial density, before switching on the speckle potential. 
A more elaborate way would be to tune the two-body interactions by using a Feshbach 
resonance [4, 37]. In this case, by scanning a magnetic field, one can continuously go from 
the independent-particle regime to the strongly interacting regime. This offers the possibility 
to study in a controlled way how localization is affected by interactions. 

In this paper we focus on the impact of disorder in the independent-particle regime 
and discard both atom-atom interactions and quantum statistical effects. Such a situation 
can be reached with ultra-cold atoms produced from a Bose-Einstein condensate or a Fermi 
degenerate gas by opening the trap and decreasing the atomic density far below the quantum 
degenerate regime. Our results therefore apply to non-interacting ultra-cold gases at low 
densities. Our paper is inspired by Vollhardt and Wolfle's seminal work on electron transport 
theory [38] and is an extension of our previously published results on localization of matter 
waves in two-dimensional speckle potentials [39]. We calculate experimentally relevant 
transport quantities, such as the scattering mean free path, the transport mean free path, 
the diffusion constant and the weak localization corrections, treating the two- and three- 
dimensional cases in parallel. Contrary to the case of ID transport studied recently by 
Sanchez-Palencia et al. [40] to which this theory can also be applied, the 2D and 3D cases 
are of special interest because 3D allows for the Anderson localization transition at a finite 
disorder strength, whereas 2D is the marginal case (the lower critical dimension) where 
analytical results are particularly important. We show that it should be possible, with the 
current experimental state of the art, to observe weak and strong locaUzation effects provided 
that phase-breaking mechanisms are under control. 

The paper is organized as follows: in section |2] we discuss the optical potential and its 
statistical properties and we derive the dimensionless Schrodinger equation which governs 
the atomic motion in the speckle potential. Sections [3] and |4] (and corresponding appendices 
containing the technical details) are intended to give an introduction to the diagrammatic 
perturbation theory. We obtain analytic expressions for the configuration-averaged propagator 
and the scattering mean free path as well as for the intensity relaxation kernel and the transport 
mean free path. Section |5] (and corresponding appendices) are devoted to the calculation of 
quantum corrections to classical transport. We derive the reduced diffusion constant in 2 and 
3D, and study the strong localization onset. 
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Figure 1. Intensity plot of a 2D speckle pattern, numerically generated as desciibed by Horak 
et al. [41]. The positions x and y are given in units of the speckle correlation length f . 



2. Atomic Hamiltonian Dynamics 

In the present section, we formulate the general description of the single-particle dynamics 
for noninteracting cold atoms in disordered speckle potentials. 

2.1. Light shifts 

When an atom is exposed to electromagnetic radiation, it is polarized, and its energy levels 
are shifted. In the case of interaction with a laser light field, these energy shifts are called 
light shifts [42]. In the dipolar approximation, the light shifts are proportional to the field 
intensity evaluated at the center-of-mass of the atom. If the field intensity is space dependent, 
so are the light shifts, and a moving atom experiences dipolar forces altering its trajectory. 
By conveniently tailoring the space and time dependence of the light field, one can produce a 
great variety of potentials for guiding the atomic motion. 

In the present paper, we consider the interaction of a two-level atom (mass to, internal 
electronic ground state |g), energy separation hut a to the excited state |e) with natural energy 
linewidth HT, electric dipole moment rfg) with a monochromatic electromagnetic laser field 
£{r) (wave number /cl, wavelength Al = 27r/fcL, angular frequency = c/cl). The two- 
level description is appropriate for atoms like strontium [43] with a nondegenerate electronic 
ground state or atoms like rubidium whose ground-state degeneracy is lifted by a strong 
magnetic field [44] . 

The coupling strength between the atom and the electromagnetic laser field is described 
by the angular Rabi frequency h^l{r) — —de£{r). The detuning (5l = wl — oja from an 
optical resonance is generally small compared to the transition frequency. Si, ^ wl- In 
this case, anti-resonant interaction terms can be ignored (rotating-wave approximation) [45]. 
Hereafter, we assume that the atom is initially prepared in its internal ground state, and that 
the laser detuning satisfies Si, ^ F, JIl, where f^L is the mean value of the Rabi frequency. 
The transition amplitude to the excited state is then small, and the ground-state light shift at 
position r caused by the laser intensity /(r) = (ieoc)|f is 



V{r) 



h\n{r)\^ nr T i{r) 



(1) 



4(5 8 (5l 
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The saturation intensity /g is a characteristic of the atom under consideration. The atomic 
motion is then governed by the effective Hamihonian H = /2m + V{r). 

The atomic evolution is purely Hamiltonian only up to some time tj: because of 
the coupling to the photon vacuum fluctuations, the light-atom interaction also contains a 
dissipative term which limits the temporal coherence of the atomic wave function. The 
average inelastic scattering rate or inverse inelastic time can be calculated using the optical 
Bloch equations [42] and reads 

7i = - ~ T — (2) 



Here Vl = ^(^) is the average light shift, given in terms of the average intensity /l = I{r): 

Coherent interference of multiple scattering amplitudes can affect the atomic dynamics if this 
dissipation is controlled at an arbitrarily low rate by using a sufficiently far-detuned speckle 
field (5l ^ r at constant average potential height Vl- 

2.2. Speckle potential 

A monochromatic speckle pattern with random intensity and phase variations defines a 
disordered potential V{r) as displayed in figure [T] The corresponding electric field 8{r) 
is a superposition of many complex field amplitudes with zero mean. Then, as stated by 
the central limit theorem, the real and imaginary parts of the field amphtude are uncorrelated 
Gaussian random variables [15]. Thanks to the Gaussian moment theorem [46], all correlation 
functions of the electric field can be expressed in terms of the pair correlation function 
'Psi't'i't'') ~ £*{r)£{r'). The ensemble average typically restores translational invariance, 
V£{r,r') ~ Vsir — r'), and one defines the dimensionless field correlation function or 
so-called complex degree of coherence [47]: 



, , £*{r + r')£{r') 
|£(r)|2 

The index o? = 2, 3 indicates that the correlation function depends on the choice of a two- 
or three-dimensional geometry as discussed below. Since the light-shift potential ([T]i is 
proportional to the square of the field, the speckle potential itself does not have a Gaussian 
distribution. Nevertheless, because of the Gaussian character of the underlying field, all 
potential correlation functions can be split down to sums of products of the field correlation 
function (|4]i. The potential pair correlation function V{r')V{r' + r) is proportional to the 
fourth-order field correlation £*{r')£{r')£*{r' + r)£{r' + r). It is given by 

V{r')V{r' + r) = V^[l + \jd{r)n (5) 
Obviously, the potential dispersion is proportional to the average potential value Vl itself. 
We rewrite the light-shift potential as V{r) = Vl[1 + (5V(r)] such that the dimensionless 
fluctuations 5V{r) have the reduced two-point correlation function 

Vd{r) - SVir')SV{r' + r) = |7d(r)p. (6) 

Vdir) is normalized such that Vd{Q) ~ 1 and decays to zero over a characteristic spatial scale 
C, known as the correlation length of the speckle potential. As shown below, is a crucial 
parameter for the dynamics of the cold atoms in the speckle potential. 
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A two-dimensional (d = 2) speckle as displayed in figure[T]can be produced by reflection 
of a laser from a rough surface or by transmission through a phase mask [48]. In the far- 
field z ^ R from the phase mask, where z denotes the distance from the phase mask to the 
observation plane and R denotes the radius of the phase mask, the speckle interference pattern 
can be regarded as quasi two-dimensional since the speckle grains are very elongated in the 
z-direction, orthogonal to the illuminated surface. Confining atoms in the plane transverse to 
z in the far-field thus realizes a 2D situation. A uniformly illuminated circular phase mask 
with radius R yields a complex degree of coherence that takes the following form for small 
relative distances r ^ z [15] 

Ji{akLr) 

72 r =2 . (7) 

Ji{x) denotes the first-order Bessel function. The geometrical factor a w R/z <C 1 is the 
numerical aperture of the imaging device; a typical order of magnitude in recent experiments 
is a w 0.1 [10] or better, a ~ 0.45 [49]. The speckle fluctuations originate from the coherent 
superposition of purely monochromatic wave vectors. This implies a diffraction limit: two 
points in the speckle field are correlated if sufficiently close to each other, limr^o 72 (r) = 1. 
In the 2D case, the correlation function 72 (r) decreases to zero over a characteristic length 
scale ( = 1/ aki, = Al /27ra such that the potential fluctuations are uncorrected only for 
distances larger than (. 

To produce a three-dimensional {d = 3) disordered configuration, the speckle grains 
should be obtained as the interference pattern of many wave vectors spanning the largest 
possible angular aperture. Ideally, this situation corresponds to the interference pattern 
obtained inside a closed optical cavity, e.g., an integrating sphere. The complex degree of 
coherence is then given by [50] 

sin(fcLr) 

73(r)^^^, (8) 

where the correlation length is now C = 1 /^l, corresponding to a numerical aperture a 1. 

Another possibly interesting configuration is provided by the speckle field at the 
proximity of a rough interface illuminated by monochromatic coherent light [51]. Due to the 
contributions of evanescent components, the average intensity then decreases exponentially 
with the distance z to the surface. If the surface is rough on short scales, the near-field speckle 
correlation length ( is smaller than Al, which is impossible for far-field speckle patterns. 
Experimentally, one would have to bring the atoms sufficiently close to the surface in a 
controllable way and restrict the dynamics to a plane parallel to the surface [52] . This situation 
deserves a special study and will not be considered in the following. 

A popular model for disordered potentials in the quantum transport literature is a 
spatially (5-correlated potential with a Gaussian distribution of potential strength. The 
previous considerations show that a monochromatic speckle potential is neither (5-correlated 
nor Gaussian. However, only the potential fluctuations as seen by the moving atoms are 
relevant for the atomic dynamics. We will see in the following that speckle fluctuations appear 
effectively (5-correlated only when the atomic de Broglie wavelength is much larger than the 
correlation length Attaining this quantum regime requires sub-recoil cooling techniques. 



2.3. Dimensionless Schrddinger Equation 



The correlation length ^ of the speckle fluctuations defines a natural physical length scale. 
In turn, it also defines natural scales for momentum, wave number, velocity, energy, angular 
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frequency, and time: 



PC 



hkc 



tq mC, 



(9) 



In the 3D-speckle case ([8]), one has C = ^ such that = hki^/m = vji is the famOiar 
atomic recoil velocity, whereas = 2Eii is twice the atomic recoil energy En = ^mv\. In 
the 2D-speckle case (|7]i, C, = {ak^)^^ and = au/^ <C f^, whereas Eq = 2a^Efi <^ E^. 
Scaling all dynamical variables with these units (fc = p/his the wave vector of the atom), 

p = r/C, r = t/ri, k = fc/fcc, e^E/E^ (10) 

the Schrodinger equation reappears in the dimensionless form idrli^) ~ H\il') with the 
HamUtonian 



T] + J]SV{p). 



In the rescaled units, the canonical commutator reads [p, k] ~ i 
strength of the potential fluctuations in units of the correlation energy 

Vl hr r /l 



E, 



C 



(11) 

The parameter t] is the 
(12) 



The constant term 77 in the Hamiltonian (fTTT i can be reabsorbed by fixing the origin of energies 
at Vl. In spatial representation and rescaled units, the stationary Schrodinger equation at 
energy e is now: 



Lo P 



v5V{p)] ^(p)-0. 



The 2-point correlation functions of the potential fluctuations in our rescaled units are 



V2{P) 



V3{p) 



p 

sin (p) 



(13) 



(14a) 



(14/7) 

P 

The equation ( fT3] l differs from the Helmholtz equation obtained for classical 
electromagnetic waves propagating in random dielectric media. In our case the potential 
fluctuations do not depend on energy. Hence the treatment for matter waves (be it atoms or 
electrons) is much simpler than for classical waves, where this energy dependence implies 
significant corrections to dynamical quantities such as the transport speed of light [25]. 



3. Effective Medium 



The states of a wave scattered by a disordered potential are different from one realization of 
disorder to another As a consequence, only expectation values obtained by averaging over 
many such realizations provide a useful characterization of the transport processes at work. 
In the language of standard quantum transport, tracing out the speckle impurities leads to a 
dispersion relation with finite spectral width 7k, for the plane-wave states \k). Equivalently, 
but put into the language of atom optics, averaging over speckle realizations introduces an 
effective medium that is characterized by a complex wave vector K{e) for the propagating 
matter wave at energy e. In this section, we calculate the dispersion relation in the weakly 
disordered regime where diagrammatic perturbation theory [20,21,53] proves particularly 
powerful. This theory is well known for the description of classical waves in a medium with a 
fluctuating index of refraction or the transport of electrons in a disordered metal. However, in 
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contrast to point-like impurities encountered for the scattering of electrons, the fluctuations of 
the optical potential exhibit spatial correlations leading to anisotropic scattering. In this sense 
the scattering of atomic matter waves in a speckle potential shares many similarities with the 
scattering of light in nematic liquid crystals where correlations between thermal fluctuations 
of the nematic director play an important role [54] . 

3.1. Retarded Green operator 

The retarded Green operator or resolvent G(e) for the stationary Schrodinger equation 
(fTSl ) at reduced energy e is the Fourier transform of the forward time evolution operator 
8(t) U{t) = 27 / de G{e) e~*^'^, with the Heaviside function Q{t). The Green operator 
satisfies the equation 

Gie) = Go{e) + Go{e)r,SVG{e). (15) 

The free retarded Green operator Go{e) = [e — Hq + iO+]^^ is diagonal in momentum 
space with propagator matrix elements {k'\ Go(£) \k) = {2ttYS{k — k')Go{k, e), where 
6{k) = S{ki) ■ ■ ■ S{Kd), and 

GoiK,e) = [e~ K'^/2 + i0+]-\ (16) 
Iteration of ( fT5] l yields the usual Born series (suppressing energy arguments for brevity) 

G = Go + rj GoSVGo + GoSVGoSVGq + ... (17) 
Taking the configuration average of the series, one obtains 

G = Go + ry^ Go 5VGa5V Go + rf Go SVGqSVGoSV Go + . . . , (18) 

where the linear term in t] vanishes because 5V ~ 0. It appears that knowing the 
average resolvent requires to calculate all higher-order correlation functions of the potential 
fluctuations SV. However, one is interested in determining the effect of the disorder on the 
energy levels £(k). To this end, the Born series is recast into the following form, known as 
the Dyson equation: 

G = Go + GoE Go + GoE GqE Go + . . . = Go + GoE G. (19) 

The retarded self-energy operator E(£) contains all irreducible correlation functions, 
i. e. correlations that cannot be split into products of independent factors by suppressing a 
single propagator Go [55]. Recognizing a geometric series, one can formally solve the Dyson 
equation as G(e) = [Go(e)~^ - E(e)]"^ 

The average over many realizations of disorder restores translational invariance. 
Consequently, G{e) is diagonal in momentum space just like the free-space propagator Go (e). 
Since the disordered potential also preserves space isotropy on average, as exemplified by the 
scalar correlation functions ^ and (O, the propagator matrix elements G(k, e) in momentum 
space can only depend on the momentum modulus k = \k\. The same conclusion holds for 
the self-energy E(k, e) such that 

G{K,e) = [e - K^/2-T.{K,e)]-\ (20) 
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3.2. Disorder-broadened dispersion relation 

Information about the relative weight, energy, and life time of excitations dressed by the 
disordered medium is contained in the spectral function [56] 

Ai.,s) ^ -2ImGK.) ^ (21) 

The spectral function is positive, A{k, e) > (because the retarded self-energy has a negative 
imaginary part), and normalized to unity J {de/2'!T)A{K, e) = 1. It can thus be seen to 
represent the probability density for excitations with wave vector k to have an energy e. Its 
trace over momentum states yields the average density of states per unit volume. 



^ A{K,e) = 27r7V(£). (22) 



The spectral function of a free particle, Aq{k, e) — 2TT5{e — k^/2), projects onto the energy 
shell = 2e, such that the free density of states reads as usual 

"^'^'^ = / (0 ~ = (iF^^e)'^/-^ (23) 

Here Sd = J dfirf = 2tt''-^^ /T{d/2) denotes the surface of the unit sphere, with the Euler 
gamma function r(a;). 

Whenever the corrections due to disorder are finite but small, A{k, e) as a function of 
e at fixed k is strongly peaked around the energy defined as a zero of the real dispersion 
relation — /2 — RcI](k, e^) = 0. A Taylor expansion to lowest order around this point, 
£ — k} j2 — RcI](k,£) = (e — £k)2'(k)^^, defines the so-called renormalization constant 
Z[n)^^ = 1— 9eRcI](K, e)!^^^ . This brings the spectral function into Lorentzian form. 

The mode |k) now represents an excitation with spectral weight Z{k) at energy e^, and with 
a finite spectral width or elastic scattering rate 

7, = -2Z(K)ImS(K,e,). (25) 

This mode describes an atom that is dressed by the speckle fluctuations and scattered into a 
different mode \k') on average after a time = l/7s- 

Equivalently, one can determine the wave vector ^(e) corresponding to a given energy 
e as the solution to the complex dispersion relation e — k^/2 — e) ~ 0. This is the 

standard approach in optics [47] and atom optics [57]. Determining the effective dispersion 
relation of matter waves in disordered speckle potentials is thus reduced to calculating the 
self-energy. 



5.5. The weak scattering approximation 

The self-energy defined by the Dyson equation ( fT9b can be expanded as a power series 

S - ^ S„, (26) 

n>2 

where each term E„ is the sum of all those irreducible diagrams that contain products of n 
field correlation functions (|4|i. 

For a weakly disordered system, the self-energy can be calculated analytically because 
only the first contribution E2 needs to be computed. Indeed, the detailed analysis in appendix 
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[Appendix A| shows that the ratio of two consecutive terms S]„ and is proportional to 

the effective scattering parameter g = rj/K. If the effective scattering is small, g ^ 1, the 
series expansion for E can be truncated after the first diagram because the effective medium 
deviates only slightly from free space. Therefore 

S«S2, (27) 
which is known as the weak scattering or Born approximation. In terms of the atomic kinetic 
energy, the weak scattering condition reads e = E / E(; ^ / 2 ^ rj^ ot 

A.!L._|_«i. (28) 

The weak scattering condition A <C 1 determines the range of validity for the diagrammatic 
perturbation theory. It can be given a simple physical interpretation by using a semiclassical 
approximation. The propagation of an incoming plane atomic wave in a speckle potential can 
in principle be computed using the Feynman path integral. If one assumes that the speckle 
potential is sufficiently weak, the so-called thin phase grating approximation [58] can be used. 
There, the classical trajectories rc{t) are unaffected by the potential (i.e. remain straight 
Unes), and the additional phase accumulated along each trajectory is simply / V{rc{t))dt/h. 
Along a path with length (, the typical accumulated phase will be V^C/hv = ^ A/2. The 
weak scattering condition A <C 1 is thus equivalent to the requirement that the accumulated 
phase be small, i. e. to the applicability condition of the thin phase grating approximation. The 
atomic wave is then only slightly distorted and scattered after traveling a distance (, which 
also implies that the scattering mean free path is larger than (. 

Equivalently, the weak scattering condition (|28] l corresponds to a small quantum 
reflection probability for a particle that is scattered by a ID potential barrier with height Vl 
and linear size C,. To see this, two cases may be distinguished: 

(1) When the potential fluctuations exceed the correlation energy (Vl > Eq or 
equivalently 77 > 1) the weak scattering condition A ^ 1 implies E ^ Vi, or equivalently 
e ^ jy; the atom flies well above the potential fluctuations. The standard quantum reflection 
coefficient [59], in our reduced notations 

4(£ - ?/) 



R 



1 



(29) 



Asin^(2V£^)_ 

then is indeed small, i? <C 1, since the oscillating term is bounded, | sin(x)/x-| < 1 with 
X = 2^e — rj. This case describes the regime of classical atomic motion. 

(2) In the opposite regime of small potential fluctuations (Vl < o'" equivalently 
77 < 1) the weak scattering condition can of course also be realized with a fast atom, again 
without quantum corrections to the classical transport. More interestingly, the weak scattering 
condition can be met even in the case E < Vl, i. e. even when the atomic energy lies below 
the average potential height. The weak scattering regime then is realized if £^ ^ i?^ where the 
large atomic de Broglie wave length AdB ^ C averages out the short-scale fluctuations and 
makes the effective disorder weak. In terms of the reflection coefficient ( [29] l, this corresponds 
to the case e < 1] and | sinh(x-) /xj w 1 with A ^ 1. We will see below that this is the regime 
where quantum corrections to classical transport become important. 

Finally, the weak scattering condition ( l28T l can be rewritten as 

E:^Ea = VllEc^. (30) 
Weak-scattering perturbation theory is valid for atoms with a sufficiently high kinetic energy 
compared to the characteristic energy i^A- We will see in section 15.41 that, in d = 3 
dimensions, E'a essentially is the mobility edge [60] which separates extended states with 
E > from localized states with E < Ej\. 
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3.4. Scattering mean free path 

Within the weak scattering approximation, the scattering mean free path, a central quantity 
characterizing the disorder, can be calculated analytically from the microscopic parameters as 
follows. The relevant self-energy contribution in the Born approximation dZTl ) is given by the 
momentum convolution of the potential correlation function Vd with the free Green function 
Go (as detailed in [Appendix A| i: 

^2{n,e)^Tf j ^^Vd{K-Ki)Go{Ki,e), (31) 

To lowest order in A = ij^/e the self-energy can be taken on-shell, E„ = Y.2{K,eK) with 
£k. = For the same reason, the renormalization constant can be approximated by 

Z{k) w 1 such that 7^ = — 2ImEK. Reverting to dimensionfull quantities, the scattering 
rate defines the elastic scattering mean free path ig/C — ^./is for a quasi-monochromatic 
wave packet centered around the momentum k — k/(: 

C 2lmT,^ 

— = -. (32) 

The scattering mean free path defines the distance over which a particle travels on average 
without being scattered. The population of the wave packet decays over the distance r by the 
factor e~'^^^% analogously to Beer's law in optics [47]. Consistently with the weak disorder 
condition jS^I ^ £, we have 7^/2 <C e and equivalently: 

k£, > 1. (33) 

Note that the scattering mean free path in the Born approximation is inversely proportional to 
A and thus to such that it does not depend on the sign of the laser detuning 

Taking the imaginary part of equation iT\\ with the help of lmGo(K, s) = — 7r5(£— k^/2) 
gives the inverse scattering mean free path ( l32b in the form 



Here, Vd{K,6) — Vdi'^n] sin||) denotes the angular correlation function as a function of 
the scattering angle 9 between k and Ki at fixed on-shell momenta k = ki = k(. The d- 
dimensional angular integration measure is dil2 — dO (integration range from to 27r) and 
dfls = 2tt sin 6 dO (integration range from to tt). 

3.4.1. 2D speckle The 2D Fourier transform of the potential fluctuation correlation ( |14a| i is 
a convolution of two identical disks: 



7'2('^) = 8[arccos|-|^l-(|) ]e(2-Ac). (35) 

In the angular correlation function 7'2(2K|sin ||), the Heaviside function 0(1 — K|sin||) 
restricts the scattering angle to |siii 1 1 < When k < 1, this condition is always fulfilled, 
and all angles are possible. When k > 1, the scattering direction is restricted to a maximum 
scattering angle \0\ < d„-^ = 2arcsin(l/K): as k grows, the differential scattering cross 
section is increasingly peaked in the forward direction. This is illustrated by figure |2(a)| a 
polar plot of the 2D phase function 
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Figure 2. 2D: (a) Phase function (36) for atomic wave vectors fc^ = 0.4, 1.0, 1.4, 1.8 (from 
left to right: blue, green, orange, red) in units of the speckle coiTelation length = l/akj^. 
Scattering is nearly isotropic for slow atoms fc^ <^ 1. For fast atoms fc^ S> 1, scattering is 
strongly peaked in the forward direction, with a maximal deviation 6m ~ . 
(b) Disorder parameter k£s as a function of the reduced matter wave number fcf for different 
disorder strengths r] = Vl/E^; = 0.2, 0.4, 0.8, 1.2, 1.6 (thin curves from left to right). The 
thick red line, connecting points of k£s where A = 1, indicates the hmit of validity of the 
weak scattering condition. 



For fast atoms k ^ 1, the differential scattering cross section is strongly peaked in the 
forward direction which clearly reveals the anisotropic nature of the scattering process. In 
this case, w 2/k <C 1. For slow atoms k ^ 1, the differential scattering cross section 
becomes isotropic. In this case, the correlation function 7^2 (p) can be approximated by a delta 
function, i. e. a constant Fourier transform 7^2 (f^) ~ 47r. Hence, the angular dependence is 
lost for K <C 1 and already the first scattering event randomizes the direction of scattering. 
This isotropic scattering limit, obtained for AdB S> C' corresponds to the s-wave scattering 
limit in scattering theory. 

Using the angular correlation 7'2(2K|sin 1 1) of ( [35] l in ( [34] i. it is possible to calculate the 
scattering mean free path. For k = 1, we can do the integral exactly and obtain 

^^^-2R^' ^C=l. (37) 

We recall that C, = {aki^)~^ is the 2D speckle correlation length and r/ = Vi^/Eq the reduced 
potential strength. It is also possible to obtain analytic results in the limiting cases k 3> 1 and 
K ^ 1 where the approximations sin a; w x and arccosa: — x\/l — w ^, respectively, can 
be made: 



kis-^jfj: fcC<l, (38a) 

k^s^"-^^, kC»l. (3Sb) 



3^(fcC) ^ 
32772 

The condition A < 1 implies the bound Ms > 1/ (27r) , such that weak scattering A ^ 1 
indeed describes weak disorder Ms ^ 1, even at very low momenta. At higher momenta 
(cf. ( |38/?t ), weak scattering A < 1 implies the bound £s > ff Ci which agrees with the 
intuitive expectation that the scattering mean free path cannot be considerably shorter than 
the 2D speckle correlation length ( itself. Figure [2(b)] shows a plot of kig as a function of fci^ 
obtained by numerical integration of ( l34l i. The boundary A = 1 indicates the limit of vahdity 
of the weak scattering approximation. 
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(a) 




(b) 



A > 1 



1.5 



Figure 3. 3D: (a) Effective phase function j40t for atomic wave vectors fcf = 
1.0, 2.0, 4.0, 6.0 (from left to right) in units of the speckle coiTelation length = I/Zcl, 
equivalent to velocity in units of recoil velocity, kC, = v/vj^. All plots with fcf < 1 are 
identical. When kC, I, the phase function is strongly anisotropic and displays the maximum 
scattering angle ~ 2/fc<^. 

(b) Disorder parameter kls as a function of the reduced matter wave number kQ for disorder 



strengths r? = VlZ-^c = 0-2, 0.4, 0.8, 1.2, 
red line connects all points of kls where A = 
scattering condition. 



1.6 (thin curves from left to right). The thick 
: 1, indicating the limit of validity of the weak 



3.4.2. 3D speckle In three dimensions, we consider a speckle pattern inside an ergodic 
cavity. In this case the 3D Fourier transform of the fluctuation correlation function (I14M 
is a convolution of two identical spherical shells such that 

T'aM = — e(2-K). (39) 

Since V3{k, 9) = Vz{'2.n sin |) diverges in the forward direction ^ 0, we plot in figure [3(a)| 
the effective phase function including the angular Jacobian, 

_ sing 7^3 

As for the 2D case, the plot shows bounded scattering 16*1 < 9^^ for fast atoms k > 1 and 
unrestricted scattering for slow atoms k < 1; exact backscattering 6 = tt appears suppressed 
due to the angular Jacobian. The inverse elastic mean free path ( l34l i for = 3 is given by 

^ - ^ [(^C)-' Q{kC - 1) + {kCr^ e(i - kO], (41) 

in terms of correlation length ^ = fc^ ^ and speckle strength 77 = Vl/S^^. 

The condition A < 1 implies the bound k£s > ^, such that weak scattering A ^ 1 
indeed describes weak disorder k£s ^ 1, even at low momenta. At high momenta, weak 
scattering A < 1 implies that k£s > f ^C' i- the lowest achievable scattering mean free 
path is of the order of the 3D speckle correlation length ( itself. Figure [3(b)| shows a plot of 
k£s as a function of k( as obtained by dTTT i. 

In summary, we have derived the elastic scattering mean free path analytically as a 
function of all relevant physical parameters. We have taken into account the correlations 
present in the speckle pattern, using a diagrammatic perturbation theory up to order in 
the reduced speckle strength. The validity of the perturbation theory is limited to the regime 
A = r/^/e = /EE^ <C 1 of not-too-large speckle fluctuations (at fixed atomic energy) and 
not-too-slow atoms (at fixed speckle strength). Not surprisingly, fast atoms are only weakly 
deviated by the disordered potential (kig ^ 1). On the other hand, we find that a strongly 
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scattering disordered medium (Mg ~* 1) can be obtained for far-detuned speckle fields of 
moderate strength and sub-recoil cooled atoms. 



4. Diffusive transport 

In the present section, we show that the matter-wave dynamics is diffusive on long time 
scales, and determine the corresponding diffusion constant by using a variant of VoUhardt 
and Wolfle's diagrammatic perturbation theory [38]. 



4.1. Density kernels and continuity equation 

In the course of its propagation in a disordered potential, the initial matter wave is rapidly 
turned into a diffuse matter wave invading the entire scattering region. The dynamics of this 
process is described by the disorder-averaged and forward-propagated probability density to 
find an atom at position p at time r > 0: 



p{p,T) = e(T)n(p,r) - e(T)Tr{e(T)n(p)} - e(T)(p|e(T)|p) (42) 

with n{p) = \p){p\ the local density operator. Propagating the initial density matrix qq 
in time, one finds p{p,t) = 6(t) (p| U{t)qqU{t)'^ \p). Its Fourier transform w) = 
J cV^p dr cxp[i((jjT ~ q ■ p)] p{p, t) is given by 



p{q,uj) = J j^-^ go{K,q)<pQ{K,q,uj), (43) 

where all information about the initial density distribution is contained in go{K,q) ~ 
(k+I go with k± = k ± q/2. The subsequent propagation is determined by the density 
propagation kernel 

<i>o(/^, g, - 1 ^ J ^ k\ q, e, Lo) (44) 
defined in terms of the intensity propagator 



^K,K\q,e,iu) = {t^'+\G{e+)\K+){K^\G^e-)\K'_). (45) 

Here, G{e+) is the full retarded propagator ( fTSl ) at energy e+ = e + uj/2, and G'^(e_) the 
full advanced propagator at energy e_ = e — uj/2. Note that the ensemble average is done 
after taking their product which means that all correlations between different amplitudes are 
included. In solid-state physics, this kernel is also known as the "particle-hole" propagator 
because both a retarded and an advanced Green's function appear. 

The average local current density is the expectation value j{p, r) = <d{T)Tr{g{T)j{p)} 
of the usual current density operator [21] j{p) — i [p«(p) + h{p)p]. Again going into 
Fourier components, it reads 

/d'^K 
j^--^go{K,q) ^iiK,,q,uj), (46) 

which is the analog of ( l44b . i.e. the convolution of the initial distribution with the current 
kernel 

^,{K,q,uj)= j ^ j ^A^'$(/^,/^',g,e,w) (47) 
that transforms indeed like a vector 
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Any Hamiltonian of the form H = /2m + V{r) leads to a so-called continuity 
equation that describes the local conservation of the probability density [21]. The only 
difference with the standard continuity equation in our case comes from the ensemble average 
that requires to impose forward-time propagation with the Heaviside distribution Q{t). The 
corresponding continuity equation 

drp{p, T) + V-j{p,T)^ S{t)p{p, 0) (48) 
features the initial condition on the right side. This conservation equation is of course valid 
for any initial distribution; in terms of the kernels ( l44l i and ( |47] i. it takes a very simple form: 

— iw(f>o('«i Q, i^) + iq ■ q, i^) = 1- (49) 

This exact relation between the intensity and current kernel does not rely on an expansion 
for small q or uj (in contrast to the approximate equation used by Vollhardt and Wolfle and 
the quantum transport literature.) This is as far as general kinematic definitions can take us. 
We now have to specify the dynamics, i.e. the equation of motion obeyed by the intensity 
propagator ( |45] ). 



4.2. Bethe-Salpeter equation, linear response and diffusion equation 

The momentum matrix elements ( l45T l define a four-point operator $(e,(jj) = 
G{e+) (g) Gt(e_) that obeys a Bethe-Salpeter equation: 

$(e, uj) = [G(e+) ® Gt (e_)] + [G(£+) ® Gt [/(e, lj) $(£, u;). (50) 

The first term on the right-hand side represents the intensity propagation in the effective 
medium with uncorrelated average propagators ( |20l i. All correlated scattering events are 
described by the irreducible intensity vertex U. The Bethe-Salpeter equation actually defines 
U, much in the same way that the Dyson equation iT% defines the self-energy E. 

Starting from this Bethe-Salpeter equation, the quantum kinetic theory described in 
Appendix B permits to calculate the current kernel (|47] | as function of the density kernel (|44] |: 



iq ■ ^i{K,,q,uj) = ^oin, q,Lu) + iT{n)q ■ k. (51) 

This expression is valid in the linear-response regime and for large distances and long times, 
qis ^ 1 and cuTg ^ 1. It features the transport time t{k) given by 



ri.)-J 2n2nme) J (2.)- ■M<s)Ai^' ,e) (1 - ■ f^'Wi^W ,s). (52) 

The spectral function A{k, e) has been defined in jTH . andJV{e) is the corresponding density 
of states. The two coupled equations ( |49] l and ( fSTT i are easily solved for the density relaxation 
kernel ( l44l i in the diffusive regime, i.e. to leading order in q and oj: 

^ol/*, q, - . X 2 - (53) 

— + lJ(Kjq'^ 
Here, the momentum-dependent diffusion constant is: 

= (54) 
Transformed back to time and position variables, the density relaxation kernel 



1 r p2 



<^o{k,p,t) ^ cxp 



(55) 



takes indeed the well-known Gaussian form that obeys the diffusion equation for a unit source 
term: 

[dr - P(k)V2]$o(k, P, t) - Sip)6{T). (56) 
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4.3. Boltzmann transport of matter waves 

The irreducible vertex operator U entering equation ( |52] | describes the average scattering of 
the local probability density. U can be expanded in a power series just like the self-energy 
(Ell): 

C/ = ^ C/n, (57) 

n>2 

Un is proportional to the speckle strength 77 to the power n, and contains all irreducible 
contributions with n field correlations and at least one correlation between the retarded 
and advanced amplitude. Due to the non-Gaussian character of intensity correlations, its 



diagrammatic representation (shown in Appendix C 1 differs from the standard form for 
Gaussian potentials. 

Generally, U cannot be calculated exactly since correlations of arbitrary order are 
involved, and one has to resort to an approximation. In the Boltzmann approximation (also 
known as the independent scattering approximation), the infinite series ( |57] | is truncated after 
the lowest-order contribution U2 such that 

C/(k, k', e) w Ub{k., k', s) = 77^ _ ^'-) (58) 



The corresponding Boltzmann intensity ( l53l ) describes multiple scattering as a sequence of 
scattering events where both retarded and advanced amplitudes travel along the same path. In 
other words, all interference effects have been discarded. This approximation thus provides 
a microscopic justification of the classical Boltzmann-Lorentz transport theory for non- 
interacting particles in the presence of quenched disorder that has been successfully applied 
to a large number of physical systems, ranging from the Drude transport theory of electrons 
in metals [26] to the radiative transfer equation in optics [27, 28]. 

The Boltzmann approximation ( fSSl l for the intensity vertex is very similar to the weak 
scattering approximation dZTl ) for the self-energy. In fact, these two approximations are 
intimately linked by the local conservation of the probability density as expressed by the 
continuity equation ( |49] l (which is guaranteed by a Ward identity ( IB. 4b in the diagrammatic 



perturbation theory, see Appendix B 1 



In the present case, inserting the vertex ( fSSl l into the general expressions (|52] | and (l54b 
defines a dimensionfull Boltzmann transport mean free path £b according to 

Vb = 1^ = -T^ = (59) 

a n a 

To zeroth order in A = v'^ /s, the spectral functions in ( |52] l reduce to the on-shell projector 

Ao{k, s) = 2TrS{e ~ such that 

where Vd{K,6) — 7'c;(2k| sin 1 1) is the angular potential correlation function with 6 the 
angle between k and k' at on-shell momenta k ^ k' ~ k(^. The transport mean-free path 
is the average distance required to completely erase the memory of the initial direction of 
propagation. It is related to the scattering mean free path ( |34] | through 



= 1 - (cos 9) = 1- dQd cos 6 UK, 6), (61) 

«B J 

where the cosine of the scattering angle is averaged over the phase function ( [36] l or ( |4Qb . This 
term is known as the anisotropy factor of the scattering process. For fully isotropic scattering, 
it is of course zero, and £b = ^s- But for strongly anisotropic scattering, {cos 9) can take a 
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Figure 4. Plot of the anisotropy factor as a function of the reduced atomic wave number k(^. 

(a) 2D = l/akj^y. For small wave numbers (fcf — > 0) the anisotropy decreases 
({cosd) 0) and the scattering becomes isotropic (£g — ► is). For lai'ge wave numbers 
(A;C ^ 1). the anisotropy becomes more pronounced ((cos 9) 1). 

(b) 3D = I/^l): Large wave numbers show strong scattering anisotropy. But even for 
small wave numbers, scattering is never really isotropic, (cos 9) = ^ for all kC^ < 1, as 
already indicated by the anisotropic phase function plotted in figure[3(a)] 



value close to 1. In this case, a large number of scattering events is necessary to deviate the 
particle completely, £b S> £s- In the language of diagrammatic quantum transport theory, 
the anisotropy factor (l6Tl i is due to vertex corrections [56] that are obtained by summing the 
entire series known as ladder diagrams; see |Appendix C| for details. 



4.3.1. 2D-speckle In the limiting cases where analytical solutions are found (cf. ( |37] | and 
(l38al l). 

fcC < 1, {62a) 

kC= I (62b) 

fcC > 1. (62c) 

A plot of the 2D anisotropy factor (cos 6) (computed numerically) as a function of kC, is 
shown in figure |43T1 For small wave numbers k( 0, one has {cos 6) and £b f-s'- 
the scattering is isotropic. When kC^ 3> 1, the ratio £b/^,s scales as (fcC)^- This can be easily 
understood because the phase function limits the angular integration to \6\ < 6'ni ^ l/^C 
such that 1 — {cos 6) ~ ^{(^'^) oc 0^. Thus roughly (fcc^)^ independent scattering events are 
needed to fully erase the memory of the initial direction. In other words, the monochromatic 
laser photons with limited projected wave vectors akh which are present in the speckle field 
can only weakly deviate atoms with large momentum k ^ aki^. 



4, 

TT^ -4 



4.3.2. 3D-speckle For the 3D case, the integration of ( l60l l yields the exact result: 

= ^ [(^C)-' e(fcc - 1) + (fcC)-' e(i - fcc)]. (63) 

In terms of the scattering mean free path (j^TJ we have 

^B = ^[(fcC)'e(fcc-i) + e(i-fcC)]- (64) 

As one can see in figure l43n a slight anisotropy I-q ~ |£s remains for all kC, < 1, implying 
(cos 9) = ^. This is explained by the absence of scattering around the backscattering direction 
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after multiplying the finite differential cross section by the vanishing angular Jacobian, as 
already evidenced by the anisotropic phase function plotted in figure |3(a)| This residual 
anisotropy is due to the long range correlations in the optical potential, as exemplified by 
the divergence of Vz{i^), see eq. (|39] l, near k = 0. At higher momenta, kC^ > 1, the ratio 
£b/^s scales as {kQ'^, for the same reason as in 2D. 

This closes our study of the classical transport properties of monochromatic matter waves 
in correlated speckle potentials. Essentially, we have found diffusive dynamics, as expected 
for particles in a conservative random potential, with a transport mean free path displaying 
a strong anisotropy for fast atoms, and becoming approximately isotropic for cold enough 
atoms. In the following, we investigate quantum corrections to classical transport. 



5. Coherent multiple scattering 

5.7. Quantum corrections to classical transport 

Within the Boltzmann approximation, all quantum interference effects are discarded. At first 
sight this seems reasonable since any such effects could be expected to be suppressed by the 
ensemble average over all possible realizations of the random potential. This means that the 
disorder average singles out products of amplitudes and conjugate amplitudes traveling along 
the same paths in the same direction where no phase differences are present. In the language of 
electronic quantum transport, these are the "particle-hole contributions". They are insensitive 
to dephasing processes and therefore correspond to classical propagation. 

It was realized however, that this argument is too simplistic [61] for phase-coherent 
systems, where interference between amplitudes of different scattering paths can occur This 
can be easily understood by considering the return probabiUty to a given point in which case 
all scattering paths are closed loops. Two waves propagating in opposite directions around 
any such loop have zero phase difference and interfere constructively (unless a magnetic 
field for charged particles is applied, or dephasing processes are at work). This constructive 
two-wave interference enhances the return probability to twice the classically expected value. 
An enhanced return probability in turn implies a reduced diffusion constant for the onward 
propagation, an effect known as weak localization. 

Quantum corrections to the Boltzmann transport picture are accounted for by including 
the sum of all counter-propagating amplitudes in the irreducible vertex U used in (|52] |: 
U = C/b + C [38, 53]. The so-called Cooperon contribution C can be expressed in terms of 
the diffusive intensity kernel ( |53] | using a time-reversal argument for the advanced amplitude 
as 

r<i I \ l-(cos6') 1 

It is a strongly peaked function around the backscattering direction Q = k + /«' = 0. 
The anisotropy factor (1 — (cos 6*)) can be justified by considering dressed Hikami boxes 
as explained in Appendix D Inserting ( |65] | into the general definition ( |52] |. one can replace 



the double integral over k' and k" by a single integral over Q (with suitable cut-offs for very 
small and very large momenta, see discussion below) such that the corrected inverse diffusion 
constant ( |54] | reads 

1 If If <^'^Q 1 \ 



V*[k,oj) Vb{h)\ nAfoie^) J (27r)'^ -icj + I?b(k)Q2 
The interference correction inside the parenthesis is the leading term in an expansion in 
powers of l/fc^e < ^/kis ^ 1. For this reason the denominator contains the free density of 
states Afo{e) instead of the disorder-averaged density of states Af{e). 
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This perturbation theory for the inverse diffusion constant 'D*{k)^^ ~ 'Db{k)~^{1 + e) 
is only meaningful for a small correction e < 1 such that the resulting diffusion constant 
is bounded from below as 'D*{k) > I?b/2. Vollhardt and Wolfe [38] have devised a self- 
consistent resummation of the perturbation series that allows in principle to reach the strong 
localization threshold 0: one simply has to replace the Boltzmann diffusion 

constant in the denominator of ( |66] | by the corrected diffusion constant 'D*{k,uj) itself. 
This prescription amounts to summing iterated loops of counterpropagating amplitudes. 
As a net result, weak localization reduces the stationary diffusion constant according to 
D* = V-Q — SV, with the quantum correction (see Appendix D for more details): 

1 f A'^Q 1 
" ^Mo{e,) J j2^W+lM' ^^^^ 
Here, the small-Q divergence of the soft mode under the integral is cut off by the real quantity 
A* := \imi^^Q[—iuj/'D*(uj)]~^/'^. The characteristic length L„ ~ ^A* can encapsulate 
several effects that limit the interference of amplitudes traveling around large loops of a 
characteristic size: (i) limited system size L, (ii) finite phase-coherence length L^, and (in) 
strong localization on a scale ^loc such that in general 

72 = 72 + 72 + 75-- (68) 

This length has been found to monitor correctly the behavior of the diffusion constant in 
bulk media (L 00) in the presence of phase-breaking mechanisms close to the strong 
localization threshold [62]. The interesting physics associated with these effects in the 
weak localization regime and at the threshold to strong localization will be discussed in the 
following subsections. 

The integral ( |67] | also diverges in the UV limit Q — ^ 00 because we took the spectral 
functions inside the double integral (|52] | at Q = 0. This divergence can thus be remedied by 
introducing an ultraviolet cutoff Qmax = C/^c given a priori by the overlap of the disorder- 
broadened spectral functions, i.e. £c = is, as discussed by van Tiggelen [31]. We have 
chosen this cut-off in our previous publication [39] . However, in the present case of correlated 
scattering, loops of counterpropagating amplitudes can only be closed on the larger length 
scale £b > is- Therefore £c = Ib seems more adequate. In particular, this choice is 
consistent with the diffusion approximation and provides particularly simple expressions in 
the following. There is no need to renormalize this microscopic cutoff length self-consistently 
as the diffusion constant, i.e. to replace £b by £*, since on the short time scale u;^^ sa tb 
during which the matter wave is scattered around a small closed loop, the weak localization 
corrections ( |66] | to classical scattering are negligible. The precise choice for the small-scale 
cutoff £c in any case can only shift the non-universal, perturbative prediction of the onset, but 
does not affect universal predictions like critical exponenents at the transition to the strongly 
localized regime; see section IST^l below. 



5.2. Weak localization 

First we will study weak localization for a monochromatic diffusing matter wave with wave 
number k (determined by the cooling technology at hand) inside a speckle field of fixed size L, 
as a function of the laser intensity /l and the detuning Sb as externally controllable parameters 
in the experiment. 

For the observation of the weak localization correction a coherent diffusive process must 
be established inside a large enough scattering region. This means that the total size of the 
scattering medium must be large enough to admit diffusion in the first place, L > ^b- And 
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Figure 5. Weak localization correction &D/ as a function of the reduced intensity /l//s 
in a speckle field of size L = 2 cm. 

(a) 2D: matter wave numbers fef e {0.8, 1.0, 1.2, 1.5, 2.0} (from left to right) with laser 
detuning (5l = 10^ T. 

(b) 3D: matter wave numbers k( e {0.6, 0.9, 1.2, 1.5, 1.8} (from left to right) with laser 
detuning <5l = 10* T. 

The weak scattering condition A < 1 is valid to the left of the thick black line. 



phase-breaking events must occur at a small enough rate 7^ = l/r^ < 7s. Otherwise, 
if the phase coherence is immediately destroyed between two consecutive scattering events, 
the propagation remains entirely classical, and the Boltzmann transport theory of section |4] 
applies. The phase coherence length = ^D-qt^ is the scale beyond which phase breaking 
mechanisms destroy weak localization. Both finite sample size and finite phase coherence are 
taken into account by using the characteristic infrared length (l68T l in the weak localization 
regime (^loc 00) as 

To ensure that interference corrections can be observed experimentally, one has to satisfy both 
the diffusive and coherent transport condition, ^ ^b- 

For atoms experiencing the light shift ([T]) inside the speckle field, the phase-breaking 
mechanism at work is inelastic photon scattering, i. e. absorption and spontaneous reemission 
of a photon into a different field mode accompanied by a recoil momentum kick for the 
scattering atom. Thus = = 7^^^, with the inelastic scattering rate 7; cx Vi^j^i, given 
by (|2l). The coherence time can be made arbitrarily large by increasing the detuning (5l at a 
fixed potential strength Vl oc /l • We will show in the following that the phase coherence 
requirements can be met with reasonable values for relevant experimental parameters. 

The 2D weak localization correction WT\ relative to the Boltzmann diffusion constant 
( |59] ) reads in dimensionfuU units and in the regime = <C given by ( 

5D__^ HL*/eB) 
Db tt Mb 



(70) 



Noticeable corrections can be expected for strong disorder where Mb is not too large. 

Figure |5j a) shows the relative weak localization correction SD/Db as a function of the 
laser intensity /l for different initial atomic velocities, at fixed laser detuning (5l = 10^ F 
and speckle size L = 2 cm. At its limit of validity A = 1, our theory predicts that the 
weak localization correction SD reaches already 20% of the Boltzmann diffusion constant 
Db itself for fcC = 2.0. For a smaller wave number fc^ = 1.5, the value of SD/Db rises to 
55%. Since C, = l/a/cL 3> I/^l, experimental evidence of 2D weak localization requires 
initial temperatures for the atomic sample well below the recoil temperature. In turn, this 
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means using a Bose-Einstein condensate as the initial atomic matter wave. As a general rule, 
the colder the atoms, the larger are the interference corrections. 

In 3D, the weak localization correction (|67| i relative to the Boltzmann diffusion constant 
( |59] | reads for > = £b 



The 3D interference corrections are small in the weak disorder regime Ms ^ 1- The 
relative weak localization correction 5D / Db as a function of the laser intensity for a speckle 
size L — 2 cm and detuning S = 10''^ F are shown in figure ISlb). As expected, the 
largest interference corrections are obtained when kC, < 1, which means initial sub-recoil 
temperatures. Here again, an experimental study may require a Bose-Einstein condensate as 
the matter wave source. 



5.3. Strong localization: 2D 

The 2D interference correction (ITQI i to the diffusion constant diverges with L* oo which 
indicates that a perfectly phase-coherent wave in an infinite disordered 2D system is in fact 
always localized, as predicted by the single-parameter scaling theory [30]. The corrected 
diffusion constant vanishes, D* 0, at the threshold determined for /c£b > k£s ^ 1 by 

HL,/£b) = ^Mb. (72) 

This condition defines a curve in ((Jlj-^l) parameter space. In figure [6(a)] we plot the 
corresponding "phase diagram" showing the boundary between the weak and the strong 
localization regime for different atomic velocities. The curves are almost straight lines 
because the Boltzmann mean free path £b in the denominator of ( fTOl i scales as rf cx V^. 
As a consequence 6D/Db scales as (/l/(5l)^, with small corrections from the logarithmic 
dependence on L^./£b- Consequently, the strong locaUzation onset has a linear dependence in 
the ((5l, -^l) plane with a slope determined by the atomic wave number k. 

For each point on these curves, one can deduce the corresponding values for the 
multiple scattering parameters. For example, for Rubidium 87 atoms {m ~ 1.44 lO^'^^kg, 
Xa = 780nm, /, = 1.67mW/cm2, F/27r = 6MHz) at kC = 1.2, L = 2cm, /l = 77/^, 
and (5l = 10^ F, we find a scattering mean free path is = 0.8 iim, a transport mean free path 
Ib =4.1 fim and a phase coherence length sa 2 mm. This places the strong localization 
threshold at k£s ~ 0.81. These numbers are of course to be taken with a grain of salt since 
they are obtained by applying the weak scattering approximation quite close to the limit of its 
validity (at the transition point A w 0.83). 

In the strong localization regime, extended atomic wavefunctions become exponentially 
localized as a function of the distance, and the corresponding localization length ^loc enters 
as a new length scale. In a bulk system L ^ oo the characteristic length ( |68] | then reads 
1/LJ = 1/L^ + l/ ^jqp [62] . Together with equation ( |72] ). this determines the 2D localization 
length ^loc as function of the atomic wave vector and the other experimentally relevant 
parameters. 

We study in figure |6(b)| the characteristic length scales on both sides of the strong 
localization threshold as a function of the laser intensity /l for kC, = 1.2, L = 2 cm, 
a = 0.1, and fixed laser detuning 5l = 10^ F. With increasing laser intensity, the phase 
coherence length (blue dashed curve) decreases since the probability of spontaneous 
photon scattering increases. The Boltzmann transport mean free path £b (turquoise dashed 
curve), a purely local quantity, shows no particular singularity, but the corrected mean-free 



21 




Figure 6. 2D: (a) Strong localization onset defined by 5D = Db in (<5l, /l) phase space for 
different atomic wavenumbers fcf € {0.8, 1.0, 1.2, 1.5, 2.0} (from right to left) in a speckle 
field of size L = 2 cm. For each value of kt^, the strong localization regime lies above and 
to the left of the coiTesponding line. The red fine corresponds to the criterion A = 1. SoHd 
curves below can reach the strong localization onset within the weak scattering regime; dotted 
curves are merely extrapolated. 

(b) Logarithmic plot of the Boltzmann transport mean free path ig, the weak localization 
corrected transport mean free path £*, the phase coherence length L^, and the locaUzation 
length 5ioc as a function of laser intensity /l in units of saturation intensity Is for a fixed 
atomic wave vector fcf = 1.2 and detuning (5l = 10'^ F inside a speckle field of L = 2 cm. 
At the strong localization thi'eshold /l = 77 /s, the connected transport mean free path 
£* = 2mD* / hk vanishes. For higher intensities the matter wave is locahzed with a finite 
localization length ^loc. 



path i* (violet dashed line) plunges to zero together with the corrected diffusion constant 
D* = ftfc^* /2m at the threshold value /l = 77 Is- This threshold value corresponds precisely 
to the transition point on the phase boundary in figure [6(a)| reached with = = 2 mm. 
At the same point, the localization length ^loc (red full curve) comes down from infinity 
and tends to the expression derived from ( |72] | in the limit ^loc ^ of strong locahzation: 
6oc = ^sexp [f fc^e]- 



5.4. Strong localization: 3D 

Contrary to the 2D case, the 3D correction ( fTTT l remains finite as oo. Therefore, a 

transition to the strongly localized regime cannot be driven solely by large-scale coherence, 
but requires strong enough local disorder This essential difference between the 2D and the 
3D case can be traced back to the fact that a random walk in 2D returns to the origin with unit 
probability, whereas in 3D it escapes with finite probability to infinity [63,64]. 

As shown by Vollhardt and Wolfle [65], the self-consistent weak localization prediction 
(6T\ allows to calculate the critical exponents of the localization length and the diffusion 
constant at the Anderson phase transition in an arbitrary dimension d > 2. We shall use their 
findings for d = 3. The integral in ( |67] | can be evaluated with an arbitrary UV-cutoff (/£c 
using the identity x^/{l + x"^) = 1 - 1/(1 + x^). At the threshold dV = Vb, it yields a 
transcendental equation for i» [21, 65]: 

= arctan -— . (73) 

Here the disorder parameter 7 and its critical value 70 are defined as 
1 
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The ratio 7/70 determines whether (|73T l admits a solution > and thus a finite localization 
length. The usual graphical solution of bx = arctan x shows that there is no finite solution 
for in the diffusive regime of small disorder 7 < 7o- In this case, the matter wave shows 
truly diffusive dynamics on all scales larger than £b with finite diffusion constant D* > 
weakly localized by the correction (ItTI ). 

In the localized regime 7 > 70, ( |73] | admits a finite solution > that tends critically 
towards infinity as 7 70. Using arctan(iH./£c) ~ 7r/2, one finds i* ^ l7^7o| 
with the critical exponent 1^ = 1. Both the diagrammatic perturbation theory and the scaling 
theory [33] thus lead to the same critical exponent 1^ = 1. 

Universal quantities like critical exponents are independent of the cut-off £c, in contrast to 
the disorder parameter and the critical threshold as defined in (l74l i. The parameter 7 simplifies 
considerably if we identify the cut-off length £c with the Boltzmann transport mean free path 
as justified in section lsTl above. With this choice 7=1/ ki-g, we locate the strong localisation 
threshold jc ~ 7o or 5D / Db = 1 at the critical disorder strength or critical loffe-Regel 
parameter [31] 



In order to estimate whether this strong localization threshold can be reached with current 
experimental techniques, we use the perturbative results of section [3] and |4] to calculate the 
laser intensity and detuning for different atomic momenta. For sufficiently large and phase- 
coherent systems, ^ £b, we have 5D/Db oc (/l/<5l)^- In the (/l, <^l) parameter plane, 
the strong localization threshold 6D/Db = 1 is thus characterized by the simple linear scaling 
II oc (5l. This can be seen in the phase diagram of the 3D strong localization onset for 
different atomic wave vectors in figure [7(a)| 

For an atom prepared at recoil momentum k = 0.9 /cl, and a realistic detuning 
(5l = 10''r, we locate the strong localization threshold at a very reasonable laser intensity 
/l — 63 /g. At this point £s = 0.09 fim, £b = 0.13 fim and the phase coherence length 
— 9 /im. Precisely at the transition point, the corrected transport mean-free path vanishes, 
£* = 0. Since ^ £b the coherent transport condition is well fulfilled. Again these numbers 
have to be taken with a grain of salt since the underlying perturbative description is bound to 
break down in the strongly disordered regime. 

The self-consistent equation ( |73] | permits to describe the locahzed regime. Close to the 
threshold, it reads 



Together with ( |68] ) this expression determines the 3D localization length ^loc. 

Figure |7(b)| shows a logarithmic plot of the relevant length scales on both sides of the 
localization threshold for a fixed atomic momentum of fc = 0.9 fcL as function of the intensity 
II at (5 = 10'' r. Just as in 2D, the phase coherence length decreases monotonically because 
spontaneous photon emission is enhanced as the power increases. The strong localization 
threshold is reached for /l = 63 Is, where the corrected transport mean free path £* vanishes 
together with the diffusion constant. As expected the localization length ^loc comes down 
from infinity at the threshold and tends toward its bulk value ^loc = |^b/(§ — ^^^b) given 
by ( |76] ) in the regime ^loc <C -L^, L. 

Often the critical behavior of is stated as a function of the energy which enters in our 
case the Boltzmann transport mean free path £b {k). For sub-recoil atoms k( < 1, the disorder 
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Figure 7. 3D: (a) Strong localization onset defined by SD = Db in (<5l, /l) phase space for 
different atomic wavenumbers fcf € {0.6, 0.9, 1.2, 1.5, 1.8} (from right to left) in a speckle 
field of size L = 2 cm. For each value of fcf , the strong localization regime lies above and 
to the left of the coiTesponding line. The red fine corresponds to the criterion A = 1. SoHd 
curves below can reach the strong localization onset within the weak scattering regime; dotted 
curves are merely extrapolated. 

(b) Logarithmic plot of the Boltzmann transport mean free path ig , the weak localization 
corrected transport mean free path £*, the phase coherence length L^, and the locahzation 
length ^loc as a function of laser intensity /l in units of saturation intensity Is for fcf = 0.9 
and S = 10* F. At the strong localization threshold /l = 63 h, the corrected transport mean 
free path £* = 3mD* /hk vanishes. For higher intensities the matter wave is localized with a 
finite localization length ^i^^ . 



parameter can be expressed in terms of the energy making use of ( l30b and (l63T l: 



where Ec denotes the mobility edge. Up to the non-universal factor 70 = \nr/3 very close to 

unity, it is given by the weak-scattering energy scale Ea defined in ( l30l l. The critical exponent 
as function of energy is still the same, cx {Ec — E)~^ with v ^1. For above-recoil atoms 
kC, > 1, the threshold condition 7 = 70 leads to an energy such that the weak-scattering 
parameter A > 1 lies formally beyond its limit of validity; in this regime, our first-order 
perturbation theory cannot give trustworthy estimates for the mobility edge. 

Similarly, we can calculate the critical exponent of the diffusion constant decreasing to 
zero on the diffusive side 7 < 70 of the transition, D* cx (70 — 7)* (x [E — EcY with the 
same critical exponent s = 1 as predicted by the scaling theory [33]. This is particularly 
remarkable since the diagrammatic perturbation theory is in principle only valid for weak 
disorder and a priori cannot describe the strong localization onset. The accurate description 
of the critical exponents within this theory however can be seen as a strong hint that the 
quantum interference corrections, which are responsible for weak locahzation, also remain 
the dominant contribution for strong localization. 

6. Conclusion and Outlook 

We have formulated a quantum transport theory for non-interacting matter waves in two- 
and three-dimensional disordered optical speckle potentials with a finite, and rather large, 
correlation length Making use of diagrammatic perturbation theory for weak scattering, 
we have determined all relevant microscopic transport parameters (scattering mean free path 
£s, transport mean free path £b and Boltzmann diffusion constant Db) that are necessary to 
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describe the average diffusion process, as functions of experimentally controllable quantities 
such as the initial atomic momentum, the laser detuning and the laser intensity. 

As a consequence of the speckle correlations we encounter an essential difference 
compared to the common treatment for Gaussian white noise. The differential single 
scattering cross section for correlated fluctuations is anisotropic and its degree of anisotropy 
depends strongly on the initial energy of the matter wave. This is a general property of 
correlated potentials which has important implications for the small perturbation parameter 
that determines the limit of the weak scattering regime. Together with the strength of the 
potential fluctuations Vi^, the speckle correlation length C, enters — in form of the correlation 
energy = h'^/mC,'^ — the important energy scale Ea = ^l/^C' which determines the 
minimum energy for weak scattering. The validity of our results is restricted to the regime 
where the atomic energy is larger than the weak scattering energy Ea- In this regime the 
atoms are only weakly scattered by the speckle fluctuations. 

Using a linear response theory, we calculate the weak localization correction to classical 
transport. These interference corrections are sensitive to phase-breaking mechanisms such as 
spontaneous photon emission. These spontaneous dissipative processes can be maintained at 
harmless rates for realistic values of the experimental parameters. We have shown that in this 
case the weak localization correction 6D can reach a considerable fraction of the Boltzmann 
diffusion constant Db within the weak scattering regime for atoms at recoil or sub-recoil 
temperatures. Even the strong localization threshold is within reach of current experimental 
techniques; we have calculated the corresponding critical disorder values and localization 
lengths. In order to measure the diffusion constant and its weak localization corrections in 
an actual experiment, one could for instance release a confined atomic cloud and monitor 
its lont-time spread inside the speckle field by time-of-flight or in-situ imaging techniques. 
The crossover to the strongly localized regime would appear as a freezing of the diffusive 
process as the cloud reaches the localization length, with a final density distribution mapping 
out the longest-lived localized eigenmodes of each realization that permit the corresponding 
localisation length to be measured. 

The extension to non-monochromatic matter wave distributions (such as the ones used 
in the recent publications [40,66]) is straightforward and will be the subject of a subsequent 
publication. 

As the use of sub-recoil temperatures favors localization effects, the most promising 
candidate for matter wave experiments seems to be a Bose-Einstein condensate. In this 
respect, another extension of this work would be to consider the dynamics of a Bose-Einstein 
condensate [40] in the 2D and 3D speckle field and to study the role of interactions (for 
instance using the approach developed in [67]) and quantum statistics on localization. Another 
possible extension of this work would be to study the influence of the internal spin structure 
of the atom. 

Finally, this work could be extended to disordered magnetic potentials. One immediate 
advantage compared to the optical speckle potential would be the Gaussian character of 
the magnetic potential fluctuations [68]. A first numerical study of the expansion of a 
Bose-Einstein condensate in a one-dimensional disordered atomic waveguide has revealed 
characteristic signatures of localization in the non-interacting regime [69]. This could open 
the door to further experiments studying disorder on atom chips [70, 71]. 
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Appendix A. Self-energy E and small perturbation parameter 



Let us recall the self-energy series 

where S„ represents all one-propagator irreducible diagrams of order n in field correlations 

and in powers of 77 = Vi^/ Eq, the speckle strength in units of correlation energy. The first 
terms are: 

E2 = •■ ■• =: ® (E) (A.2) 

E3 = . (A.3) 

E4 = (g)- — (g)— — (g) — -(g) + g)- — g)- ®—® 

+ • . . . + •—•^—^ + • . . . (A.4) 



E5 = (g— •— •- — -• -g) + + . . . (A.5) 

In these diagrams, every straight line represents the free propagator Gq defined in ( fTSI l. 
"Irreducible" means that by cutting any such single propagator line, the diagrams do not split 
into independent parts. A dotted line like • • = 7/7^ represents the field correlation function 
(|4|. A double dotted line • • = 'rf'Vd = |7d P represents the intensity correlation function 
(|6l), that we note rf'Pd = ® ® for simplicity in the following. To evaluate a diagram, one 
writes down its expression in terms of correlation functions and propagators and integrates 
over all free internal variables. This integration is most conveniently done in momentum 
space since the free propagator Go{k, e) is then diagonal. 

As an example, let us calculate the first diagram E2 = ® — «> explicitly. The potential 
fluctuation operator g) = rjSV is diagonal in real space and therefore translation invariant 
in momentum space: {k'\SV\k) = SV{k — k'). The Fourier transform of the potential 
fluctuations is defined by 

6V{k) = J d^p 5V{p) e-'^ P, SV{p) = J ^ SVin) e'^ f. (A.6) 

The 2-point correlation function of the potential fluctuations in momentum space reads 

5V{ki)5V{k2) = {2^Y5[ki + k^) Vdif^i), (A.7) 
where Vdi^) is the Fourier transform of Vd{p) that was defined in (|6]l. The matrix elements 
of E2(£) are given by {k'\ E2(e) \k) = {2tiY5{k — k') E2(/t, e); as expected, the self-energy 
operator is diagonal in momentum space and isotropic. Its diagonal entries are 

E2(K,e) = ry2 /" ^^Vd{i^- Ki)Go{n^,e), (A.8) 
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which is simply the momentum convolution of the potential correlation function Vd with the 
free Green function Gq. 

The appearance of odd terms S2f/+i in the formal series jA.lj reflects the non-Gaussian 
character of the potential fluctuations. We now wish to determine the effective small parameter 
and corresponding range of validity of the self-energy series. To this aim, we discuss 
separately the high-energy and low-energy cases n = \/2e '3> 1 and k ^ 1, respectively. 

In the high-energy limit, it is possible to determine the dependence of I]„ on the 
ultraviolet momentum k 1 by simple power counting. Each irreducible diagram 
contributing to I]„ contains rt — 1 internal propagators Go(e), p field-correlation functions 
{0 < p < n), and (n—p)/2 intensity-correlation functions. Taking into account all momentum 
conservation laws, this leaves exactly (n + p)/2 independent variables Ki that can be chosen 
to be the arguments ^d{i^i) or Vdii^i) of the correlation functions ( I14ab and (I14/?b . Because 
of the strict momentum cutoff, these correlation functions constrain the norms Hi to remain of 
order unity or smaller The only dependence on n comes from the Green functions Go ( Hm , £) 
that are evaluated at momenta k,„ ~ n J^i'^i'^i with coefficients S {0,±1} that 
describe the diagram's topology. Linearizing around the on-shell value e = k^/2 for Ki, 
each Green's function contributes a power k^^ such that I]„ ~ a„7]"K^^" = ann{r] / k)"^ , 
where a„ is related to the number of n-point irreducible diagrams. We thus can identify the 
effective expansion parameter g = that should be small, which in turn justifies the choice 
of the weak scattering condition A = rj^ /e = 2g^ ^ 1 in section [33] above. Of course, 
such an analysis of a superficial degree of divergence [72] can only determine how the small 
parameter depends on energy, but not fix numerical factors of order unity. We can rewrite the 
ratio of the self-energy jA.lb to the kinetic energy e as 

^ — ' K 

n>2 

In this high-energy regime, non-Gaussian terms n = 2q + 1 contribute to the series. In other 
words, fast atoms live in a potential that has definitely non-Gaussian statistics beyond the 
leading-order Born approximation n = 2. 

The low-energy limit k <C 1 in 2D allows for a similar analysis. Each of the {n + p)/2 
correlation functions in S„ tends towards its finite limit as k ^ and thus becomes effectively 
^-correlated in real space. The contributions from the (n + p)/2 integration measures 
and from the n — 1 Green's functions yield a scaling for S„. If n is even, 

Tjn is therefore dominated by the diagrams with p ~ 0, containing only intensity-correlation 
functions and diverging like If n is odd, the dominant contribution comes from the 

diagrams with the smallest number of field correlations, p ~ 3, that appear first in S3 
and reappear subsequently in higher non-Gaussian terms 'S2q+i- One can then rewrite the 
dominant contributions to (lA.lb in the form 

S/e - ^{a2qg^'^ + a2q+iK^g^'^+^), (2D and k < 1). (A.IO) 

Remarkably, g = tj/k is still the expansion parameter, and remains small if ij ^ k ^ 1. 
Terms which are neghgible compared to both their neighbors in the series can be omitted. This 
applies to all odd terms when K^g^'^^^ ^ 5^*^^^, i- e- g at fixed g. In this low-energy 

and weak-scattering regime, the self-energy does no longer depend on the field-correlation 
functions, and only the Gaussian terms survive: 

S/e - ^02,5^'', (2DandK3<5). (A.ll) 

9>1 
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The 3D case requires a separate discussion. As far as the pure intensity correlations 
(p = 0) are concerned, actually the same reasoning as for 2D holds because the 1ow-k- 
divergence of the intensity correlation oc is compensated by a supplementary factor from 
the integration measure. However, the pure field-correlation functions 73(k) = S{1 — k) 
cannot tend to a constant as k ^ 0, but project all integration momenta onto the unit sphere. 
Consequently, a small-K-contribution of the integration measure and the correlation functions 
can only come from the (n — p)/2 variables intervening in the intensity-correlation functions: 
it reads k"^^. The contribution of a propagator is either (as in the 2D case) if it does not 
depend on a field-correlation momentum, or independent of k in the limit k ^ because field 
correlation momenta of order unity remain present. Thus, the diagrams depending only on the 
field-correlations {n = p) behave like k'-\ At fixed n and p, the dominant diagrams mixing 
field- and intensity-correlation functions are those with the largest number of propagator lines 
independent of the p field-correlation variables. This happens when field correlations can 
be written as products of the largest number of independent field-correlation sub-diagrams 
that never cross the intensity-correlation lines. An example for such a diagram is the first 
contribution to S5 shown in ( IA.5I ) that displays a Sa-type field correlation inside a S2 intensity 
correlation. In all dominant cases, the field sub-diagrams will contain at most 3, 4 or 5 
vertices (as those shown in E3, S4, and E5), because higher-order field correlations could 
be factorized into these elementary ones, thus yielding an additional independent propagator. 
Writing p ~ 3ni + 4^2 + 5713, where the n,; are non-negative integers, the largest possible 
number of sub-diagrams is obtained by maximizing the sum ni + + n^. The number of 
propagator lines giving a contribution to the diagram is then n — p + ni + n-z + — 1, 
and the total contribution of the diagram is ,{2-n+ni+2n2+3n3 turns out that when n is 
even, the dominant contribution k^^" to I]„ comes from the intensity-correlation diagrams 
(p = 0). When n is odd, the main contribution «;^~" is due to the diagrams with ni = 1, 
i.e.p = 3. Similarly to the 2D case, the dominant contribution to the self-energy (lA.lb is 



As in the 2D case, g ~ tj/k is the expansion parameter. Again, terms which are negligible 
compared to both their neighbors can be omitted. This applies to the odd terms when 
^ ?7 ^ K ^ 1 at fixed g, i. e. in the quantum regime discussed at the end of section [33] 
In this regime, an effective (5-correlated Gaussian potential is recovered. 

Since the number of diagrams a„ grows factorially with ri, we face the well-known 
troublesome fact that even for g <C 1, our global weak scattering condition covering both the 
high- and low-energy regime, the series (lA.lb formally diverges. It can only be understood 
as an asymptotic series that can be accurately approximated by just the first few terms [73]. 
When the effective coupling constant g is sufficiently small, then a truncation to the first term 
already gives a good approximation to the self-energy. The weak scattering condition is thus 
given by g ^ 1, or equivalently hy rj k, which amounts to E ^ Ea = V^^/E^. 

Self-energy diagrams like — «i and ®— «> — 8 — «i and similarly nested higher-order 
diagrams with an outer correlation function can be self-consistently summed up to give 
the diagram ®— ®, where the thick line represents the average Green function. The same 
procedure can be applied to all diagrams like 181^ — 0—0 — -g) with two outer correlation 
functions yielding the self-consistent diagram «>— ®— «i, etc. 




(A.12) 



9>1 
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Appendix B. Diffusion from quantum kinetic equation 

This appendix is devoted to the derivation of the diffusion constant using a quantum linear 
response theory that we believe to be a slightly improved version of VoUhardt and Wolfle's 
approach [38]. 



Appendix B.l. Quantum kinetic equation and continuity equation 

We first rewrite the Bethe-Salpeter equation ( fSOl l for the intensity propagator <&(/«, /«', q, £, lo) 
in a more useful way. For the scalar quantities the product of the average propagators 
G'^{eS) ® G{e+) can be reformulated in momentum space by using the identity G* G = 

(G^ -G)/(G'^ -G^^^) as 

^ -AG(<g s a.) ^^ ^^ 

Lu — q ■ k' — AT,{k' ,q,e,uj) 

where k'^ = k' ± q/2 and AG{K,q,e,uj) ~ G'(K+,e+) — G*(K_,e_) as well as 
AY,{k, q,e,uj) = £+) — S* £_). Multiplying the denominator to the other side 

leads to the quantum kinetic equation [2 1 ] 

[uj-q- k'] k\ q, e, uj) ^ -{2'kY5{k - k')AG{k' , q, e, lu) + C[$]. (B.2) 

This is the standard form of a kinetic equation with a Fourier-transformed drift derivative 
on the left hand side, and on the right hand side first a source term and then all scattering 
information contained in the linear collision functional 



-AG{k', q, e, u) J K",q, e, uj)U{k\ k' , q, e, u). (B.3) 



C[$] = AE(k', q, e, uj)^{k, k' , q, e, iv) 

As shown by Vollhardt and Wolfle [38], the irreducible vertex U and the self-energy S 
are intimately linked through the Ward identity [21, 38] 

AS(/^',q,e,c^)- j ^^AG{K\q,e,uj)U{K",K',q,e,u). (B.4) 

Notably, the proof works also for correlated potentials and thus applies without restriction to 
our case (in contrast to scattering by resonant potentials where corrections to the Ward identity 
renormalize the transport speed of the propagating wave [25]). This Ward identity plays the 
role of an optical theorem; everything that disappears from the forward propagating mode is 
scattered into other modes. With its help, the scattering functional can be rewritten as 

C[$]= j [AG{K\q,e,uj)<^{K,K' ,q,e,uj) 

-AG{k\ q, £, ijj)^{k, K.",q, e, lu)] U{k", k' , q, e, lu). (B.5) 

By parity, it vanishes identically under integration over k', J d'^K'Cm = 0. Moreover, one 
can verify that / deAG{K, q, e,uj) = 27ri by going back to the definition of the time-evolution 
propagator and observing that U{t = 0) = 1 while defining 0(0) — |. By integrating the 
quantum kinetic equation over k' and e, we thus recover the continuity equation ( |49l ) exactly 
to all orders in q and lu. This is slightly different from Vollhardt and Wolfle's approach, 
followed by the entire quantum transport literature, which uses different kernels defined by 
sums over k' and k (instead of k' and e) that only obey an approximate continuity equation. 
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Appendix B.2. Linear response 

Starting with the continuity equation, we now need a second equation that expresses the 
current $i [n, q, oj) as function of a small gradient in density, thus defining a linear-response 
coefficient. Multiplying the QKE ( IB. 21 ) hy q ■ k' and integrating over e and k' gives an 
equation for the current kernel $i: 

ujq-^i{K,q,uj) - J ^ j ^j^{q- K.'f<^{K,K',q,e,uj) 

This is not a closed equation for $o ™d $i because it involves higher-order moments of the 
full distribution If one continues to write equations for the new unknowns, one generates a 
hierarchy that must be truncated judiciously in order to extract the density kernel $o- 

In the long time and long-distance limit q,uj ^ 0, we can make two assumptions that 
facilitate the derivation of the transport time in the linear response regime. Firstly, we perform 
the usual expansion 



l + ±i^'.q)iq.n") 



<P{K,K",q,e,uj) (B.7) 
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into the lowest-order terms in powers of the angular dependence on q ■ k', isotropic and 
vectorial. The factorized spectral function A{K',e) encapsulates the dependence on the 
modulus of k' which is constrained to the disorder-broadened energy shell; indeed, the 
Bethe-Salpeter equation ( |50l ) shows that the intensity relaxation kernel $ is proportional 
to AG{k' ,q,e,uj) which is itself very well approximated by the strongly peaked spectral 
function A{K,e) since AG(/«,e) = 2ilmG{K, e) = —iA{K,e). Secondly, by symmetry 
in K and k', the intensity relaxation kernel $ is also proportional to A{K,e). As a sharply 
peaked function of e, the kernel can therefore be well approximated by k" , q, e, lo) k, 
A{k, e) J ||-<i>(/«, k", q, e', w). Using this in ( IB.7I ) leads to the ansatz 



~ , A{K,e)A{K',e) 
^{k,k ,q,e,uj) = 



^o{K.,q,uj) + ^{k' ■ q)q ■ ^i{K.,q,uj) 
2e 



(B.8) 
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This expression in terms of the relevant kernels $o and $i is custom-tailored such that the 
consistency relations $o = *&o and = $i hold. We can therefore insert $ w $ into 
iB.6i and calculate the two remaining terms. On the left-hand side, only the $o-term gives 
a contribution since the ^>i-term vanishes by parity. On the right-hand side, the collision 
integral of the isotropic $o-term vanishes by parity if AS, AG, and U are evaluated at q = 
and uj = 0. This is allowed to leading order in q, lo because the collision integral is akeady of 
order l/r^. The remaining integral over the $i-term defines the inverse transport time ( |52] i. 
Rearranging the equation, we can solve for g • $i and obtain (ISTI i. 

Together with the continuity equation (|49] l this provides a closed set of equations for $o 
and $1 with the following solutions for ujt{k) <C 1: 

1 — iT(hi)q ■ K , ^ 'D(K)q^ 

<Poi^,q,u;)= . ^ 7 ^ , iq-^,{K,q,u)^ . ^ ^ (B.9) 

with the reduced diffusion constant = K^T{K)/d. Inserting these solutions into the 

expansion dB.Sb yields the corresponding diffusive intensity propagation kernel 

2ttJ\1 (e) —lu + V{K,)q^ 
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Appendix C. Intensity diagrams and Boltzmann transport theory 



The Bethe-Salpeter equation (pOb can be formally recast into the form 

$ = [Gt G] + [Gt ® G]i?[Gt ® G], (C.l) 

where the reducible vertex is i? = U + U [Gt (g) G] R. The first terms of the power series 
U = X^„>2 Un for the irreducible intensity vertex are the following diagrams; 

«) 

U2 = : (C.2) 



Us = ■ , + '\ (C.3) 



U4= X + « \ + +... (C.4) 







As before, "irreducible" means that these diagrams do not fall apart into independent sub- 
diagrams by cutting one of the thick lines that represent the average propagators G. The upper 
lines belong to retarded amplitudes G while the lower entries belong to advanced amplitudes 
G^^. Dotted lines connecting two (g) (resp. •) represent an intensity (resp. field) correlation 
function. In addition to the familiar potential correlations as in U2 and U4, we find field- 
correlation diagrams as in U3 and all higher orders, since the potential fluctuations do not 
obey Gaussian statistics. 

Through the Ward identity ( IB.4l i. every self-consistent diagram in the expansion of the 
self-energy corresponds to a set of diagrams in the expansion of U. For example, the diagram 
(IC.2b is linked to the self-consistent diagram g)— igi, exactly as in the electronic case. This 
means in particular that the Boltzmann approximation U ~ Ub = U2 and the weak scattering 
approximation S w S2 are consistent, and the Ward identity (IB.4l i appears as the expression 
([34I 1 of the scattering mean free path in terms of the Boltzmann vertex (ISSl l. The next two 
diagrams ( IC.3b are linked to the field correlations defining S3, and the following 

(IC.4b are linked to (»— ®— ®— ®, etc. 

In the Boltzmann approximation U ^ U2, the reducible vertex at the core of the intensity 
propagator jC.lb becomes _Rb = C/2 + U2[G'' ^ G]U2 + ■ ■ ■ = U2 + L. Here L denotes the 
so-called dijfuson, the letter L referring to the ladder topology of diagrams: 



L 



(C.5) 



We will need an analytical expression for the diffuson in the diffusive regime. To this end, we 
consider the geometric series for <I>b in the Boltzmann approximation: 



$B = + + (C.6) 



The first term on the right-hand side describes ballistic propagation, whereas the second 
term describes single scattering in the effective medium. These two terms can be neglected 
compared to the remaining diffuson contribution to the series which is given by dB.lOl ) 
with T) = V-Q. Making use of the identity dB.lb in the long time and large distance 
limit, i. e. for a; = and q = 0, the spectral functions in (IB. 10b can be expressed as 
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A{K,e) = G* {k, e)G{K, e)/Ts. Both spectral functions ^(k, e) andA{K',e) each contribute 
a factor I/tj. The products of the retarded and the advanced Green functions as functions of 
K and k' can be identified with the entering and exiting propagators to the diffuson dC.Sb in 
(IC.6b . The diffuson is therefore given by 

, 1 1 — iT{K)q ■ {k + k') 
L(k,k ,q,e,uj)= . : — . . — . (C.7) 

Only the prefactor differs form the diffusive intensity propagation kernel dB.lOl l. In particular, 
the diffuson displays the same diffusive pole as the diffusive intensity propagation kernel. 

Appendix D. Weak localization correction for anisotropic scattering 

The weak localization correction to the transport time ( l52b is contained in the intensity 
scattering vertex U = Ub + C, where C incorporates the full cooperon contribution 

C = Ca + Cb+Cc (D.l) 
The momentum diagrams corresponding to these three contributions are 



C= + + =: (D-2) 



Special attention has to be paid to the cooperon contribution in the case of anisotropic 
scattering. The simple substitution U = Ub + Ca which is sometimes used in the literature 
(cf. [53]) is not sufficient in this case. For electron transport in highly anisotropic systems 
Wolfle and Bhatt [74] have performed the calculation of the cooperon contribution for the 
conductivity tensor. In this calculation the two additional conductivity diagrams Cb and Cc 
corresponding to dressed Hikami boxes [20] are taken into account. This calculation for the 
anisotropic conductivity tensor is similar to the calculation of the weak localization correction 
to the diffusion constant that we wish to perform. We therefore need to include the same 
additional Hikami diagrams Cb and Cc adapted to our case. 

The total transport time including the weak localization correction takes the form: 

tb tb , tb f de A{k, e) f d'^n' d'^K" 



'(k) tX{k) t| J 27r27r7V(e) J {2^^ {2tt)'^ 

X [[G(At')]'G^('t') [G(k")]'G^(k") + G(k')[G^('«')]' GiK")[G^{K")f 

d'^Q 



x{l + k'- k") Vd{K' - k") J ^-^ Lo{k, Q, u;) (D.3) 



(27r)" 

where the energy dependence of e of the Green functions has been left out for better visibility 
and where 

In order to evaluate the integral over k,' and k" in ( ID. 31 ) one makes use of the peaked structure 
of the Green functions as a function of k' and k" around the point = \/2e for weak 
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disorder. To leading order in the small disorder parameter 1 / eTj <JC 1 the integral over k! and 
ft" then yields: 



TB 

t*{k) 



= 1 + 



27r27r7V(£) J {2ttY 



Loin, Q,u!) 



2[/ 



2,2l 



The function /" "(r^) is defined by [20]: 

J (Ts j = 27rA/o « 



(n + TO-2)! 



(D.5) 



(D.6) 



(to - l)!(n- 1)! 

Reinserting the functions /"'™(ts) into dP.St . the term in brackets yields the factor 
47rA/'oTg (1 — (cosO)). Evaluating the integral over e to leading order in I/eTs, we obtain 
the final result 



TB 



1 



d^Q 



1 



t*{k) ^ ' TTj\fo{e^)J (2^)^ -ilu + Vb{k)Q^ 
that leads to the corrected inverse diffusion constant ( 



(D.7) 
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